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1. 


INTRODUCTION 


This  report  addresses  the  problem  of  utilizing  spatial-hyperspectral  imaging 
capabilities  of  space-borne  sensors  to  detect  and  characterize  regions  of  atmospheric 
turbulence  and  cirrus  cloud  clutter  which  may  impact  employment  and/or  performance  of 
space  based  laser  and  airborne  high  energy  laser  weapons  systems. 

Our  concept  for  worldwide  detection,  characterization  and  mapping  of  atmospheric 
turbulence  and  cirrus  clouds  involves  use  of  satellite-borne  (and  possibly  airborne)  spectral 
and  hyperspectral  imagers  operated  in  the  UV  to  MWIR  spectral  range  in  virtual  triangulation 
geometry.  Spectral  and  hyperspectral  imagery  allows  altitude  sounding  of  atmospheric  clutter 
from  turbulence  and  cirrus  clouds.  Triangulation  geometry  allows  precise  altitude  selection  by 
cross  correlation  of  the  backscatter  signals.  The  combination  of  altitude  and  Fourier-space 
background  spectral  discrimination  will  provide  an  altitude  resolved  measurement  of 
atmospheric  clutter  from  clear  air  turbulence  and  from  cirrus  clouds,  both  of  which  may  affect 
performance  of  the  SBL  (Space  Based  Laser)  and  the  Airborne  Laser  (ABL)  systems. 

Figure  1  illustrates  the  virtual  triangulation  concept  by  forward-looking  and  backward¬ 
looking  imaging  sensors  from  a  satellite. 


Figure  1  Conceptual  Optical  Triangulation  Geometry  to  Isolate  Atmospheric  Turbulence,  Characterize 
Ground-based  Clutter,  and  Detect  Weak  Targets,  The  virtual  triangulation  geometry’  allows  mapping  of 
random  atmospheric  clutter  and  CAT  vs.  altitude  by  cross  correlation  of  forward-viewing  and  backward¬ 
viewing  imaging  signals  and  assists  in  separating  weak  airborne  targets  from  the  atmospheric  background 
clutter. 
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We  anticipate  that  the  conceptual  satellite  sensor  can  be  tested  on  ground-based  and  airborne 
platforms  prior  to  finalization  of  a  satellite  sensor  design.  This  phase  one  SBIR  project  has 
initiated  such  conceptual  tests  but  further  work  is  required,  especially  with  respect  to  airborne 
tests  of  the  concept  versus  clear  air  turbulence  (CAT). 

Anticipated  commercial  applications  of  this  research  are  automated  background  scene  content 
identification  and  weak  target  recognition  by  small,  real-time  programmable  space-borne 
sensors  and  remote  sensing  of  atmospheric  structure  which  leads  to  a  capability  for  mapping 
CAT  (clear  air  turbulence)  from  airborne  or  space-borne  detectors. 

2.  STATEMENT  OF  PROBLEM  AND  MEASUREMENT  CONCEPT 

2.1.1  Technical  Problem  Summary 

The  technical  problem  addressed  by  this  SBIR  program  is  to  devise  unique  spatial-spectral 
satellite-based  observation  techniques  to  detect  weak  targets  in  the  presence  of  background 
clutter.  The  weak  targets  that  we  address  are  in  order  of  observability: 

•  small,  low  observable  airborne  vehicles, 

•  thin  atmospheric  cirrus  clouds, 

•  atmospheric  clear  air  turbulence,  and 

•  turbulent  structure  in  the  troposphere,  stratosphere,  and  lower  thermosphere. 

This  report  concentrates  on  the  more  difficult  “targets”,  atmospheric  turbulence  and  clear  air 
turbulence  (CAT)  layers.  The  rationale  for  measurement  of  the  atmospheric  quantities  is  that 
this  class  of  backgrounds  forms  the  limiting  cases  of  background  clutter,  which  if  measurable 
by  a  satellite  passive  sensor  ensures  that  other  airborne  targets,  will  be  detectable. 

2.1.2  Technical  Approach 

We  validated  the  virtual  triangulation  -  spectral  -  spatial  weak  target  detection  and 
characterization  technique  by  analytical  simulation  techniques  and  by  experimental  simulation 
using  ground-based  spectral-spatial  structure  sensors  and  by  using  the  GLO  sensors  carried  on 
the  Space  Shuttle1. 
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Validation  bv  Simulation 


We  used  the  NSS  (Non-stationary  Stochastic  Structure)  Model2  and  LAMSS3  (Low  Altitude 
Mesoscale  Stochastic  Structure)  models  to  simulate  structure  in  the  mesopause  region  (~90 
km)  and  in  the  troposphere  (<  20  km)  respectively.  These  simulations  were  generated  using 
an  geometrical  driver  for  the  sensor  LOS  to  generate  an  artificial  data  stream  appropriate  to 
the  experimental  validation  geometries. 

Validation  bv  Experiment 

Two  experiments  were  conducted  to  validate  the  satellite  sensor  cross-beam  concept:  a  low 
altitude  experiment  measured  turbulent  structure  and  winds  in  the  boundary  layer  between 
about  50  and  300  feet  altitude.  We  intended  that  the  daytime  measurements  would  be 
coordinated  with  meteorological  tower  measurements  of  winds  and  turbulence  but  this 
experiment  could  not  be  consummated  under  the  resources  of  the  project.  A  high  altitude  set 
of  measurements  was  conducted  at  nighttime  used  the  airglow  emission  in  the  80  to  130  km 
altitude  range  to  provide  the  signal.  Daytime  experiments  used  the  Kestrel  AIRCAM  sensors, 
operating  in  the  450  to  950  nm  spectral  range.  The  nighttime  high  altitude  structure 
experiments  used  the  Fourier  Transform  Hyperspectral  Imager  FTHSI  and  supporting 
imaging  instrumentation  from  LANL  and  University  of  Arizona. 

2.1.3  Results  Summary 

The  Phase  I  simulations  and  ground-based  experiments  validated  to  a  certain  extent  the  ability 
of  a  satellite-borne  passive  atmospheric  structure  sensor  to  provide  worldwide  measurements 
of  atmospheric  turbulence  and  wind  structure  especially  thin  cirrus  clouds  and  clear  air 
turbulence.  Validation  experiment  results  may  be  applied  to  further  experiments  to  determine 
atmospheric  structure  phenomenology  and  clutter  maps  for  the  Airborne  Laser  System  (ABL), 
the  Space-based  Laser  System  (SBL),  CAT  maps  for  the  civilian  sector,  and  detection  weak 
targets  through  efficient  mitigation  of  of  atmospheric  background  structure-induced  clutter. 
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3. 


CONCEPT  VALIDATION 


The  atmospheric  structure  sounding  concept  was  validated  by  two  techniques.  Simulation  of 
atmospheric  structure  and  synthesis  of  atmospheric  backscatter  signals  for  low  altitude 
structure  and  atmospheric  airglow  signals  for  high  altitude  structure  was  carried  out  for  a 
limited  number  of  cases  and  scenarios.  Experiments  designed  to  simulate  possible  satellite 
and/or  airborne  sensor  environments  and  their  respective  signals  generated  by  structure  were 
conducted  to  investigate  the  utility  of  Hyperspectral  Imagery.  Data  from  the  GLO  sensor  that 
has  been  repeatedly  flown  on  the  Space  Shuttle  also  contributed  to  our  experimental  database. 

3.1  VALIDATION  BY  SIMULATION 

3.1.1  Simulation  of  Winds  and  Turbulence 

Atmospheric  turbulence  was  simulated  in  two  dimensions  by  the  technique  demonstrated  by 
Strugala  and  Sears,  et.  al.  (See  reference  ).  Two  different  two-dimensional  spectra  were 
adopted:  a  k“5/3  power  law  approximates  the  turbulence  spectrum  for  an  optically  thick 
medium,  and  a  k  power  law  spectrum  approximates  an  optically  thin  turbulent  medium. 
Both  spectra  are  of  the  modified  Kolmogorov  type  with  a  correlation  length  of  10  cm,  about 
the  value  expected  for  the  lower  atmosphere. 

Equation  1  illustrates  the  functional  form  of  the  spectra  for  a  two  dimensional  scene: 

P(kx,  ky)  =  [T(  p  )  T(  p )  L«  Ley  ]/  {kx/LCx} 2  +{ky/Lcy)  2]}p  (1) 

Where, 

p  =  5/6  or  8/6  for  k'5/3  or  k'm  spectral  forms  respectively, 
n  =  (p-l)/p, 

ky,y  is  the  spatial  frequency  in  either  x  or  y  dimension 

f(  n)  is  the  gamma  function  with  argument  determined  by  the  spectral  index 
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The  background  spatial  scenes  were  implemented  using  the  method  of  Strugala  and  Sears 
wherein  a  random  Gaussian  number  field  is  created,  then  convolved  with  the  desired  2-D 
spectral  function  (equation  1)  to  achieve  a  correlated  random  background  intensity  field. 
Figure  2  illustrates  the  two  dimensional  scenes  generated  by  the  simulation  routine.  Here,  we 
used  the  same  random  number  set  in  two  dimensions  to  generate  the  arrays  and  vary  only  the 
power  spectral  density  slope.  A  64x64  pixel  sample  was  cut  from  each  256x256  scene  to 
show  that  the  higher  slope  (k-8/3)  has  much  more  high  frequency  structure  as  indicated  by  the 
density  of  the  contour  lines.  Each  picture  was  digitized  to  a  six-level  contour  interval. 


Figure  2  Contour  maps  of  2-D  scenes  generated  from  the  same  Gaussian  random 
number  field,  but  with  spectral  indices  of  k  -5/3  (left)  and  \C*3  (right)  which  physically 
corresponds  to  Kolmogorov  turbulence  in  an  optical  thick  and  optically  thin  medium. 

In  order  to  simulate  the  effects  of  winds  on  the  turbulence  observations,  we  cut  subimages 
from  the  256x256  images  and  displaced  them  spatially  by  a  constant  number  of  pixels  in  the  x 
and  y  dimensions.  The  analysis  algorithms  then  were  applied  to  compute  the  complex  cross 
power  spectral  densities  of  the  spectral  images. 


The  basic  algorithms  used  for  the  triangulation  geometry  as  well  as  the  successive  frame 
analysis  in  the  vertical  direction  is  the  cross  power  spectral  density  algorithm  group  first 
applied  to  geophysical  phenomena  by  Gossard  et.  al.4  The  Cross  Power  Spectra  for  either  a  set 
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of  one  dimensional  vector  quantities  VI,  V2,  or  a  sequence  of  two  dimensional  images  II,  12 
may  be  generically  expressed  as: 

XP  (VI  i(,  V2  i,  or  II  y,  12  y)  =  CFT(V1  +iV2,  or  I1HI2)  (2) 

Where  the  data  vectors  or  image  fields  have  indices  i,  and  i,j  respectively 
CFT  is  the  complex  Fourier  Transform  with  the  vector  or  image  quantities  combined 
as  a  complex  argument. 

Whereas  Gossard  et  al  described  the  process  in  terms  of  the  cross  correlation  functions  of  a 
data  series,  modem  algorithms  such  as  the  Fast  Fourier  Transform  allow  direct  computation 
of  the  data  set.  This  analytical  procedure  was  implemented  by  using  standard  complex  Fourier 
Transform  Algorithms  contained  in  Mathcad.5 

3.1.2  Simulation  of  Triangulation  Geometry 

We  simulated  the  triangulation  geometry  in  two  dimensions  corresponding  to  the  right  angled 
intersection  of  the  low  altitude  structure  and  wind  experiment  described  later. 

3.1.3  Signal  Analysis  of  Simulations 

We  analyzed  the  simulation  experiments  for  two  cases.  The  2-D  single  beam,  vertical  looking 
geometry  that  corresponds  to  the  low  altitude  experiments  measures  MWIR  scattered  sunlight 
as  a  signal  as  the  turbulent  scattering  irregularities  are  convected  through  the  field  of  view  of 
the  sensors.  The  2-D  low  altitude  triangulation  geometry  corresponds  to  the  scattered  sunlight 
signal  in  the  region  of  intersection  of  the  sensor  fields  of  view. 

In  the  case  of  the  vertical  viewing  sensor  simulation,  we  computed  the  complex  cross  spectral 
information  in  both  real  and  imaginary  regime  and  using  standard  signal  processing 
techniques  (see  reference  4)  derived  the  AS  (Amplitude  spectra),  PSD  (power  spectral 
densities),  IPSD  (imaginary  components  of  PSD),  and  the  PHSD  (Phase  Spectra)  for  each 
pair  of  images.  Then  we  looked  for  peaks  in  the  phase  spectra  (PHSD)  which  might 
correspond  to  the  motion  of  the  turbulent  scattered  scenes  across  the  sensor  field  of  view. 
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The  two  dimensional  PSD’s  and  Phase  Spectra  (PHSD)  for  an  image  array  A  are  defined  by 
the  standard  expressions  involving  the  periodogram  PA  of  the  array  A  which  is  derived  from 
the  “cfft”  complex  Fourier  Transform  of  the  array  A.  In  the  simulation  case,  the  image 
arrays  A  are  the  computed  Gaussian  matrices  GM. 


PA  =  cfft(A),  PSD  =  |PA|2  (3) 

PHSD  =  -  arctan  (Im(PA)/(Re(PA)  (4) 


As  a  diagnostic  we  compare  in  Figure  3  the  ideal  Kolmogorov  k"  PSD  input  to  the 
program,  the  computed  periodogram  for  the  convolution  of  the  GM,  Gaussian  correlated 
noise  matrix,  and  the  raw  periodogram  of  the  array  GM  itself.  The  periodograms  will  be 
converted  to  power  spectral  densities  (PSD’s)  later  in  the  analysis 


Figure  3  Left  Image  Compares  the  Correlated  Gaussian  Noise  Matrix  1-Dimensional 
Periodograms  and  the  input  data  (smooth  line).  SME  is  the  input  PSD,  PDSE  is  the 
Periodogram  of  the  Noised  Gaussian  Matrix,  and  PDGM  is  the  Periodogram  of  the 
Gaussian  Noise  Matrix  before  convolution  with  SME.  Values  of  PDSE  and  SME  are 
normalized  to  their  maximum  values.  The  right  image  shows  the  1-D  PSD’s  (lines)  after 
a  Parzen  filter  is  applied.  The  data  points  correspond  to  the  noisey  periodogram  in  the 
left  image.  The  x  axes  correspond  to  spatial  frequency  in  units  of  inverse  pixels. 


This  exercise  confirms  that  the  PSD  of  the  simulated  random  turbulence  field  reproduces  the 
input  spectrum. 


Also  it  is  necessary  to  validate  that  phase  spectra  can  be  produced,  using  the  Gaussian 
random  structure  simulation.  Figure  4  illustrates  the  phase  spectrum  over  the  range  of  zero  to 


10 


1/64  pixel  spatial  frequency  range.  Histograms  indicate  that  the  phase  spectral  field  is  a 
correlated  Gaussian  distribution  as  expected. 


Figure  4  Phase  Spectrum  for  Gaussian  Correlated  Simulation  Scene.  The  spectrum 
extends  from  zero  to  1/64  pixel  spatial  frequency  in  both  spatial  dimensions  and  from  - 
1  to  1  radian  in  phase. 

The  Cross  Power  Spectral  Density  and  Cross-Phase  spectra  are  used  to  estimate  the  spatial 
displacement  of  the  images  and  the  vector  velocity  of  the  passage  through  the  simulated  two 
fields  of  view,  from  successive  analyses  of  the  sub-image  pairs  that  are  separated  by  20 
pixels. 

The  2-D  Cross  PSD  is  smoothed  by  the  Parzen  window  technique  and  then  converted  to  a  2- 
D  correlation  function  my  means  of  a  complex  FFT.  Figures  5A  and  5B  illustrate  the  2-D 
PSD  as  a  2-D  plot,  and  as  a  pair  of  1-D  cuts  along  the  axes. 


1 1 


SXPA 


Figure  5A  Smoothed  Cross  PSD  as  Figure  5B  Smoothed  Cross  PSD 

a  surface  plot.  Note  the  occurrence  shows  spectral  peaks  on  the  frequency 

of  several  spectral  peaks  x-y  diagonal  as  anticipated.  The  other 

peaks  appear  to  be  extraneous  features. 


The  coordinates  of  Figure  5A  are  in  units  of  spatial  frequency  divided  by  the  pixel  dimension 
of  the  image  (ie  128  pixels  with  2  pixels  being  the  Nyquist  frequency).  Units  of  5B  are  in 
inverse  pixels  such  that  0.5  is  equivalent  to  the  Nyquist  frequency  of  1/(2  pixels). 


After  conversion  of  the  2-D  Cross  PSD’s  to  cross  correlation  plots,  it  is  necessary  to 
threshold  the  data  to  bring  out  the  pertinent  features  on  the  desired  scale.  Figures  6A  and  6B 
show  respectively  the  un-thresholded  2-D  cross  correlation  plot  and  the  thresholded  image. 
The  threshold  chosen  for  this  case  was  7  times  the  standard  deviation  of  the  cross  correlation 
scene  after  removal  of  the  random  energy  component  at  zero  delay. 

The  axes  are  in  units  if  1/pixels,  hence  the  value  0.5  is  equivalent  to  2  pixel  displacement  and 
0. 1  is  equivalent  to  10  pixels  displacement.  The  length  of  the  thresholded  correlation  value  in 
Figure  6B  shows  that  the  width  of  the  correlation  results  equals  about  20  pixels  and  is  in  the 
diagonal  dimension  as  specified  by  the  original  image  pairs. 
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Figure  6A  2D  Cross  Correlation  Plot  Figure  6B  2-D  Correlation  Plot 

Unthresholded  Thresholded  at  7  Sigma 


Figures  7A  and  7B  show  the  phase  spectrum  thresholded  at  two  levels,  guided  by  the 
histogram  of  phase  (not  shown  here).  Spatial  frequencies  are  in  units  of  inverse  pixels. 


0  0.004  0.008  0.01  0.016 


TAPHS 

Figure  7A  Exceedance  Map  for  1  sigma 
Thresholded  Phase  Difference  Spectrum. 


TAPHS 

Figure  7B  Exceedance  Map  for  3  sigma 
Thresholded  Phase  Difference  Spectrum. 
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VALIDATION  BY  EXPERIMENT 

3.1.4  Low  Altitude  Wind  and  Turbulence  Experiments 

Low  altitude  wind  and  turbulence  experiments  were  conducted  at  Kestrel  Corporation  using  a 
set  of  filtered  CCD  cameras.  These  cameras  are  normally  operated  as  an  airborne  remote 
sensing  system  by  Kestrel  Corporation  and  provide  digital  images  with  8  bit  resolution  and  a 
252  x  252  image  pixel  frame  size.  The  goal  of  these  experiments  was  to  measure  solar 
backscatter  from  aerosols  and  Rayleigh  scatter  in  the  lower  troposphere,  in  the  near  vicinity  (a 
few  hundred  meters)  of  the  sensors.  Appendix  1  summarizes  the  sky-scatter  structure 
experiments  conducted  by  Kestrel  Corporation.  In  summary,  we  obtained  individual  green 
and  blue  images  with  unpolarized,  horizontally,  and  vertically  polarized  images,  plus  dark 
response  on  each  focal  plane  (cover  on),  and  uniformity  of  the  focal  plane  responsivity  by  use 
of  a  sunlit  diffusion  source. 

Figures  8  A  to  8C  illustrate  luminance  images  (combined  green  and  blue  filters)  used  to 
determine  the  focal  plane  background  dark  noise  (8  A),  the  focal  plane  response  to  a  diffuse 
illumination  source  (8B),  and  the  response  to  the  sunlit  sky.  In  the  case  of  the  sunlit  sky  figure 
8C  illustrates  the  solar  scatter  symmetry  very  well  in  both  green  and  blue.  Unfortunately  black 
and  white  reproduction  is  unable  to  show  these  features  clearly. 


Figure  8A  Dark  Focal  Plane  Figure  8B  Diffusely  Figure  8C  Luminance  image 

Illuminated  Image  of  sky  showing  green  and 

blue  scatter  from  the  sun. 

Figures  9A  and  9B  show  the  solar  scatter  pattern  observed  by  camera  V2  observed  in  green 
and  blue  channels  respectively.  Image  analysis  can  be  accomplished  using  the  composite  G+B 
to  enhance  signal  to  noise  ratio  or  the  separated  G  and  B  channels. 


Figure  9 A  Solar  Scatter  Observed  in  Figure  9B  Solar  Scatter  Observed 

Green  Band,  Sensor  V2  Blue  Band,  Horizontally  Polarized 

Several  diagnostic  configurations  were  explored  for  this  experiment.  The  dark  current 
statistical  noise  characteristics  for  the  cameras  were  measured  by  covering  the  instruments 
and  recording  data.  The  focal  plane  structure  noise  was  measured  by  putting  an  (almost) 
uniform  diffusive  surface  over  the  lenses  and  recording  data.  Finally  measurements  were 
made  in  several  wavebands,  and  in  two  polarization  modes  for  vertically  viewing  sensors  and 
for  the  virtual  triangulation  geometry  with  the  sensors  separated  by  66.5  feet  and  the  altitude 
of  intersection  about  100  feet.  Table  1  summarizes  the  focal  plane  and  objective-space 
characteristics  of  the  visible  and  NIR  cameras.  For  more  details  of  this  experiment  refer  to  the 
Kestrel  data  summary  report  in  Appendix  1 . 


Table  1 

Camera  Characteristics  and  Projected  Objective  Space  Quantities 


Camera 

Band 

Pixel  IFOV 
mrad 

FOV 

mrad 

Pixel  Footprint 
At  100  ft 

VSSHHj 

isssawMi 

V2 

Green31ue 

(Polarized) 

.91 

442  x 1075 

~  2.8  cm 

768x486 

V4 

Green,Biue 

(Polarized) 

.91 

442  x 1075 

~  2.8  cm 

768x486 

m 

MWIR 

3-5  mm 

1.4 

680  x  680 

~  4.7  cm 

486  x  486 

In  addition  to  the  visible  G-B  data  sets  an  Amber  near  infrared  (3-5  |im)  camera  was 
operated  to  view  zenith  in  the  vicinity  of  the  visible  beam  intersection.  This  camera  was 
used  primarily  to  see  if  the  signal  to  noise  ratio  on  a  cooled  MWIR  system  could  detect 
structure  in  the  lower  atmosphere.  Because  of  the  8-bit  digitization  level  (i.e.  256  intensity 
levels)  we  believed  that  this  sensor  might  be  marginal  for  detecting  atmospheric  structure  but 
the  presence  of  water  vapor  and  nearly  subvisual  clouds  may  have  enhanced  our  ability  to 
detect  turbulence  and  winds.  Figures  10A  -  IOC  compares  one  data  frame  obtained  by  the  IR 
camera  with  the  focal  plane  response  to  both  cold  and  hot  plates  placed  over  the  objective 
lens.  It  is  apparent  that  a  strong  asymmetry  in  focal  plane  response  exists. 


Figure  10A  IR  Data  Image  Figure  10B  IR  Hot  Plate  Image  Figure  IOC  IR  Cold 
Plate 

The  2-dimensional  spatial  PSD’s  (Power  Spectral  Density  estimates)  of  truncated  images 
(128x128)  pixels  and  2-d  correlation  functions  were  computed.  The  images  were  truncated  to 
ensure  that  radial  variation  in  focal  plane  noise  and  responsivity  could  be  minimized.  This 
procedure  also  (nearly)  eliminates  the  requirement  for  flat-fielding  the  data  frames. 

The  PSD’s  computed  for  the  Kestrel  AIRCAM  imager  in  measurement  mode,  dark  current 
mode,  and  diffuse  scatter  mode  are  compared  in  Figure  1 1A  and  1  IB.  Figure  1 1 A  compares 
the  measured  sky  background  spectrum  with  example  spectral  slopes,  k''  ‘ ,  k'  and  k  .  The 
experiment  spectrum  drops  into  the  noise  level  at  a  spatial  wavelength  about  30  pixels. 

Figure  1  IB  shows  that  the  focal  plane  pattern  and  dark  current  noise  (symmetric  in  k-space) 
and  the  readout  noise  (1/f  or  (k'2)  spectrum  in  the  readout  direction  may  be  a  significant 
deterrent  to  measurement  of  very  small-scale  structure  inasmuch  as  the  structure  spatial 
frequency  may  lie  beyond  the  noise  cutoff  frequency. 
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AIRCAM  DATA  PSD’S 


log(s) 

Spatial  Frequency  -  (s/128  cy/pix) 


AIRCAM  CALIBRATION  DATA 


Figure  11A  AIRCAM  sky  data  1-D  PSD 
compared  to  slopes,  -5/3,  -2,  -8/3. 

Note  that  this  camera  was  uncooled. 


IR  CAMERA  -  DATA 


log(s) 

Spatial  Frequency  -  (s/128  cy/pix) 


Figure  12A  IR  Camera  Data  PSD’s 
Compared  to  model  PSD  slopes. 


Figure  11B  AIRCAM  Calibration  PSD’s 
for  dark  counts  and  uniform  illumination 
The  readout  noise  is  the  k-2  slope 


IR  CAMERA  CALIBRATION 


log(  s) 

Spatial  Frequency  -  (s/128  cy  /pix) 


Figure  12B  IR  Camera  PSD’s  on 
Calibration  Plates  compared  to  models 


Figures  1 1  and  12  clearly  show  the  role  of  focal  plane  noise  and  readout  noise  in  the 
AIRCAM  and  IR  sensors.  The  readout  noise  is  clearly  1/f  (in  amplitude)  and  is 
approximately  the  same  magnitude  for  the  IR  sensor  sky  data  as  for  the  hot  plate  and  cold 
plate  calibration  frames.  In  fact,  the  mean  value  of  sky  background  in  the  IR  band  was  lower 
than  either  the  hot  or  cold  calibration  frame  values,  indicating  that  the  sky  was  colder  than 
either  calibration  temperature.  For  these  reasons  we  do  not  attempt  to  correct  the  IR  camera 
data  for  noise  and  cannot  pursue  analysis  in  terms  of  structure  motion  and  correlation  length. 

We  pursued  estimation  of  the  phase  drift  and  the  correlation  function  parameters  for  the 
AIRCAM  data  based  upon  the  nearly  2  order  of  magnitude  signal  to  noise  ratio  (  accounting 
for  focal  plane  plus  readout  noise  spectra).  Figure  13  illustrates  the  average  autocorrelation 
function  for  the  five  frames  of  data.  Figure  14  presents  the  “auto-Phase  spectrum  for  the 
same  set  of  AIRCAM  data  frames. 


TRP 


Figure  13  Autocorrelation  Function  Figure  14  Average  Phase  Spectrum  for 

for  5  successive  AIRCAM  Image  Frames  Same  set  of  AIRCAM  Images 

The  autocorrelation  data  clearly  shows  a  symmetric  pattern  characteristic  of  the  Gaussian 
(presumably)  characteristic  of  atmospheric  low  altitude  turbulence  as  averaged  over  5 
successive  image  frames.  The  phase  data  are  ambiguous,  but  seem  to  indicate  definite  phase 
peaks  at  selected  spatial  frequencies  which  may  be  emphasized  by  our  thresholding  method. 
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The  Parzen-  windowed  power  spectra  illustrated  in  Figure  15A  are  used  to  compute  the  cross 
correlation  functions  for  successive  image  frame  delay  times  as  illustrated  in  Figure  15B. 


CROSS  PSD  -  2  FRAME  DELAY 
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Figure  15A  Cross  PSD  for  Successive 
AIRCAM  IMAGE  Frames.  Straight 
Lines  are  K-5/3,  k-2,  and  k-8/3  slopes. 


( center  ( TRXCF4) ) 

Figure  15B  Cross  Correlation  Function 
For  Same  Image  Frames.  Displacement 
From  center  indicates  a  2  pixel  lag  in 
Diagonal  direction  (=2  mrad). 


The  results  shown  in  Figure  15A  and  B  show  that  valid  cross  PSD  and  correlation  functions 
can  be  computed  for  atmospheric  turbulence  on  a  small  scale  despite  the  high  noise  levels  in 
the  sensor  and  backgrounds.  The  2  pixel  offset  in  the  correlation  function  center  corresponds 
to  a  horizontal  “wind”  displacement  of  the  turbulent  structure  of  about  40  cm  corresponding 
to  a  wind  velocity  of  about  5  to  10  mph.  This  is  consistent  with  local  observations. 


3.1.5  High  Altitude  Experiments 

Two  types  of  high  altitude  simulation  experiments  were  conducted.  We  analyzed  the  GLO 
spectrometric  data  from  the  University  of  Arizona  instruments  flown  on  Shuttle  Flight  STS- 
85.  The  GLO  instruments  were  flown  in  the  look-forward,  look-back  configuration  with  a  + 
65  degree  azimuth  with  respect  to  the  perpendicular  to  the  spacecraft  velocity  vector  to 
provide  virtual  triangulation  geometry.  The  instrument  lines  of  sight  intersected  the  limb  at 
about  58  km. 
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The  second  experiment  was  observation  of  the  airglow  emissions  in  the  visible  to  Near-IR 
spectral  range  using  the  Kestrel  Fourier  Transform  Hyperspectral  Imager  (FTHSI).  In  both 
experiments,  the  atmospheric  airglow  emission  was  the  target  signal.  The  goals  of  these 
experiments  were  first,  to  determine  if  the  virtual  triangulation  geometry  actually  could  be 
implemented  by  space-borne  detectors  and  second,  to  determine  if  an  FTHSI  instrument 
could  produce  useable  atmospheric  data  with  sufficient  spatial  and  spectral  resolution  to 
identify  well  known  atmospheric  airglow  features. 


3.1.5.1  GLO  Airglow  Data  Analysis 

Dr.  Lyle  Broadfoot  of  University  of  Arizona,  Lunar  and  Planetary  Institute  has  conducted 
GLO  experiments  on  several  STS  flights.  A  selection  of  data  from  STS  -  85  were  provided  to 
this  project  under  subcontract.  The  geometry  of  the  STS  -  85  experiment  is  similar  to  that 
shown  in  figure  1  except  that  the  two  sensors,  GLO  5  and  GLO  6  are  oriented  to  intersect  the 
limb  with  a  line-of-sight  tangent  altitude  of  59  km.  Appendix  2  contains  selected  details  of 
the  STS  -  85  GLO  experiments.  Table  2  summarizes  the  spacecraft  parameters  and  camera 
details  relevant  to  this  experiment 


Table  2 

Spacecraft  Parameters  and  Camera  Characteristics 


Spacecraft: 

STS  -  85 

Altitude: 

400  km 

Orientation: 

Pitch  =  180  degrees,  Yaw  =  270 

Velocity 

-  7  km/sec 

Hvperspectral  Cameras: 

GLO  5  &  GLO  6 

Camera  Spectral  Range: 
Camera  Orientation: 

Elevation  Angle  (depression)  22  degrees,  azimuth  =  +  65 

degrees 

LOS  Tangent  Altitude  (center) 

59  km  (drifts  with  time) 

Range  to  Tangent: 

2120  km 

LOS  Velocity  at  tangent 

~  5  km/sec 

Geometrical  Time  Lag  of  5-6  LOS  -420  sec.  (=  2120  km/5  km/sec  ) 

Spatial  footprint 

5  km  at  earthlimb 

Figure  16  summarizes  the  GLO  airglow  data  for  one  satellite  pass.  Figures  17-21  show  the 
detailed  GLO  measurements  of  four  airglow  species  emissions  where  the  time  variability 
between  GLO  5  (forward-looking)  and  GLO  6  (backward-looking)  and  the  time  delay 
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between  the  data  sets  is  established  by  matching  intensity  features.  Each  plot  shows  a  delay 
from  the  forward  to  backward-pointing  sensors  appropriate  to  the  airglow  enhancements  at 
their  emission  altitude.  The  thin  airglow  emitting  layers  are  intercepted  by  the  sensors 
operating  at  59  km  tangent  altitude  at  two  positions,  a  near  field  between  the  tangent  point 
and  the  sensors  and  a  far  field  intercept  point  on  the  far  side  of  the  59  km  tangent  point.  The 
emission  “onion  skin”  characteristic  leads  to  two  possible  lag  times,  corresponding  to 
features  in  the  near  and  far  fields.  It  is  unlikely  that  the  sensor  vertical  resolution  is  sufficient 
to  spatially  resolve  the  far-field  component  so  we  will  concentrate  on  near  field  results. 


We  computed  the  delay  times  for  a  number  of  emission  altitudes:  85,  95,  130  and  240  km  to 
compare  with  the  GL05-6  delay  observations.  Table  3  summarizes  the  airglow 
characteristics  and  near-  and  far-field  ranges  for  the  sensor  to  intercept  the  emitting  layer. 
The  corresponding  delay  times  compared  with  the  GLO  values.  We  adopted  an  LOS  velocity 
of  6  km/sec. 


Table  3 

Airglow  Delay  Predictions  and  Observations 


Emitter 

Species 

Band 

nm 

Altitude 

km 

Range 
km  -Near 

Range 
km  -  Far 

Delay 
min.  -near 

Delay 

min.-far 

GLO 

Observation- 

min. 

— 

OH 

762 

85 

670 

6100 

1.9 

17 

7.2 

Na 

589.3 

95 

470 

6300 

1.3 

18 

4.8 

02 

762 

95 

470 

6300 

1.3 

18 

2.8 

OI 

557.7 

-130 

340 

6400 

1. 

18 

2.8 

OI 

630 

-240 

320 

7000 

0.9 

19 

2 

Comparison  of  the  computed  estimates  of  delay  time  vs.  GLO  measurements  shows  that  most 
of  the  observed  delay  times  are  about  a  factor  of  2  greater  than  the  predicted  estimates.  Some 
of  the  factors  creating  this  discrepancy  may  be  due  to  airglow  horizontal  structure,  and  the 
exact  altitude  of  the  layer  as  the  two  sensors  scan  the  same  region.  There  also  is  a  trend  in  the 
GLO  observations  to  much  higher  values  of  delay  time  in  the  OH  and  Na  emissions.  Both  of 
these  emissions  are  highly  structured  spatially  and  temporally  so  may  not  indeed  be 
statistically  stationary  over  the  period  of  the  measurement,  about  30  minutes. 
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3e0535  0H(5,1) 


j - GLO  e_58fcm - GlQ5_-5acm  [ 


3e053&  OHC5J1 


OLO  5  S?X7T~  ■  —  -  GLO 


Figure  17  OH  (5-1)  Meinel  Band  Intensity  vs.  time  from  GLO-5  and  GLO-6.  The  top 
panel  shows  the  data  at  the  time  of  acquisition.  The  bottom  panel  shows  the  same  data 
with  the  trailing  sensor  plot  time  offset  for  a  best  fit  of  the  intensity  features 
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3e0535  _Q2(0J)) 


3e0535  02(0,0) 


j - GIO  5_59km - GLO  &-4  8  min  | 


Figure  18  02  Atmospheric  Band  Intensity  vs.  time  from  GLO-5  and  GLO-6.  The  top 
panel  shows  the  data  at  the  time  of  acquisition.  The  bottom  panel  shows  the  same  data 
with  the  trailing  sensor  plot  time  offset  for  a  best  fit  of  the  intensity  features 
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3e0535  Na(5892) 


3e0535  Na(5892) 


Figure  19  Na  589  &  589.6  nm  Intensity  vs.  time  from  GLO-5  and  GLO-6.  The  top  panel 
shows  the  data  at  the  time  of  acquisition.  The  bottom  panel  shows  the  same  data  with 
the  trailing  sensor  plot  time  offset  for  a  best  fit  of  the  intensity  features 
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3e0535  01(5577) 


3e0535  01(5577) 


Figure  20  01  (557.7  nra)  line  Intensity  vs.  time  from  GLO-5  and  GLO-6.  The  top  panel 
shows  the  data  at  the  time  of  acquisition.  The  bottom  panel  shows  the  same  data  with 
the  trailing  sensor  plot  time  offset  for  a  best  fit  of  the  intensity  features 
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3e0535  01(6300) 


3e0535  01(6300) 


Figure  21  OI  (630  nm)  line  Intensity  vs.  time  from  GLO-5  and  GLO-6.  The  top  panel 
shows  the  data  at  the  time  of  acquisition.  The  bottom  panel  shows  the  same  data  with 
the  trailing  sensor  plot  time  offset  for  a  best  fit  of  the  intensity7  features 
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3. 1.5.2  Hyperspectral  Imager  Analysis  of  Airglow  Spectra 


The  Kestrel  Fourier  Transform  Hyperspectral  Imager  (FTHSI)  was  used  to  estimate  the 
capabilities  of  this  instrument  to  measure  high  altitude  airglow  spectra  in  the  400  to  900  nra 
range  from  the  ground.  This  instrument  was  not  designed  to  perform  this  measurement 
mission,  but  we  decided  to  attempt  the  measurement  and  to  test  the  focal  plane  cooling 
requirement  required  to  detect  airglow. 

Figure  22  illustrates  the  night  airglow  data  obtained  by  the  FTHSI  operating  in  an  out-of¬ 
design  dry-ice  cooling  mode.  The  cross  track  data  were  integrated  to  obtain  a  higher  signal  to 
noise  ratio  to  suppliment  the  cooling  effect.  We  found  that  the  streetlight  and  other  urban 
signals  such  as  Na  and  Hg  lines  were  observed,  but  also  we  could  identify  most  of  the 
emission  features  as  airglow  features.  Figure  22  also  identifies  the  airglow  features  including 
Na  and  Hg  emissions  from  streetlights.  Only  a  few  radiance  enhancements  were  not 
identified,  but  this  is  common  in  an  urban  background  environment.  The  integration  time  for 
these  measurements  was  order  1  second  or  less,  (down  to  TV  30  Hz  sample  rate). 

Figure  23  shows  night  airglow  features  measured  many  years  ago  as  published  in  Chamberlin6 

We  conclude  that  the  FTHSI  and  similar  imaging  spectrographic  instruments  have  the 
sensitivity  to  detect  airglow  emissions  in  the  visible-near  IR  bands  and  that  useful  space-borne 
measurements  can  be  obtained  provided  that  the  focal  planes  are  cooled  sufficiently. 
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wavelength 


Figure  22  Airglow  Spectrum  Measured  by  FTHSI  instrument  The  emission  features  are 
identified  by  comparison  with  streetlight  (Na  and  Hg)  emissions  which  also  serve  as 
wavelength  calibrations  and  by  comparison  to  well  known  airglow  features. 
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4. 


SENSOR  AND  SPACECRAFT  CONCEPT  DEVELOPMENT 


In  order  to  provide  an  initial  feasibility  check  regarding  incorporating  this  instrument  concept 
into  a  small  satellite  bus,  we  initiated  a  simple  satellite  design  concept  based  upon  the 
dimensions  of  the  MightySat  II  bus.  The  purpose  of  this  design  effort  was  to  ascertain 
whether  off  the  shelf  CCD  cameras  and  hyperspectral  imagers  could  be  coupled  with  a 
MightySat  II  class  buss  and  preserve  the  telescope  resolution,  instrument  placement,  and  other 
parameters  required  for  the  virtual  triangulation,  hyperspectral  measurement  requirements. 
Figures  24  and  25  illustrate  the  design  concept  drawings. 


Figure  24  Overall  line  drawing  of  conceptual  sensor  and  satellite  bus  based  upon 
MightySat-II  bus  dimensions.  The  instrument  dimensions  are  semi-scale,  based  upon 
commonly  available  cooled  CCD  focal  plane  detectors  for  UV  and  visible  spectral  range. 
An  MWIR  instrument  may  be  somewhat  larger  and  will  require  a  second  beam  splitter. 
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Panel  PI  shows  an  overall  view  of  the  MightySat  II  bus,  solar  panels,  and  attached  optical 
system.  P2  shows  orthogonal  views  of  the  concept  satellite  MightySat  II  bus.  P3  shows  the  third 
orthogonal  view  with  the  solar  panels  partly  folded.  The  optical  system  is  designed  to  fit  within 
the  envelope  of  the  folded  solar  panels,  with  the  exception  of  the  heat  exchanger  on  the  space¬ 
viewing  side  of  the  satellite.  P4  shows  the  overall  envelope  of  the  proposed  optical  system  the 
configuration  for  which  is  shown  in  P5.  We  plan  to  use  dielectric  layered  beamsplitters  to 
separate  UV,  visible  and  IR  optical  signals. 


Figure  25  shows  the  optical  layout  of  a  two  camera  (UV  and  IR)  sensor.  The  dimensions  of  the 
components  are  approximately  to  the  scale  of  commercially  available  visible  and  UV  cameras. 
A  cooled  MWIR  camera  such  as  those  made  by  AMBER  would  be  considerably  larger,  but  the 
height  of  the  optical  module  would  fit  would  be  scaled  to  accommodate  the  specific  unit 
selected. 


Figure  25  Telescope  mirror  layout  for  conceptual  instrument  design.  The  instrument  uses  an 
off-axis  telescope  in  order  to  maximize  cross  track  field  of  view. 
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We  also  assembled  a  very  semi-scale  model  of  a  satellite  to  show  that  commonly  available 
cooled  CCD  detectors  would  fit  on  a  MightySat-II  bus.  Figures  26A  and  26B  show  a  rendition 
of  the  satellite  concept  as  a  third-scale  model. 


Figure  26A  Front  view  of  concept  satellite  model  based  upon  MightySat  II  dimensions.  The 
optical  sensor  unit  is  a  small  module  between  the  telescope  and  the  satellite  bus. 


Figure  26B  Rear  view  of  concept  satellite  model  as  above.  The  focal  plane  and  telescope¬ 
cooling  grid  is  attached  to  the  space-directed  black  surface  of  the  MightySat-II  bus.  The 
present  concept  involves  using  heat  pipes  to  cool  the  telescope  and  focal  planes. 
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5. 


COMMERCIAL  APPLICATIONS 


Commercial  and  military  applications  for  the  multispectral,  virtual  triangulation  remote 

sensing  technique  lie  in  the  following  areas: 

•  Airborne  sensing  of  atmospheric  clear  air  turbulence  (CAT)  and  small  scale  turbulence  as 
applied  to  airborne  and  space-borne  laser  weapons  systems, 

•  Satellite-borne  mapping  of  atmospheric  turbulence,  waves,  and  winds  as  applied  to 
transmission  of  energy  by  means  of  waves  and  turbulence  from  the  troposphere  through 
the  stratosphere  into  the  mesosphere, 

•  Satellite-borne  sensing  of  atmospheric  turbulence  regions  which  would  affect  airborne  and 
space-borne  laser  weapons  systems. 

Our  commercialization  plan  incorporates  the  following  elements  which  will  be 
implemented  in  the  Phase  II  program: 

•  Modification  of  the  existing  space-qualified  GLO  sensor  system  developed  by 

University  of  Arizona,  by  improving  the  field  of  view,  pixel  resolution,  and  inclusion 
of  a  Fourier  Transform  Hyperspectral  imager  in  place  of  the  dispersive  imaging 
spectrometer.  We  also  will  incorporate  the  quad  C-40  processor.  Both  the  FTHSI  and 
the  Quad  C-40  will  be  tested  on  the  MightySai  n.  1  satellite  testbed.  This  instrument 
will  operate  in  the  UV  through  NIR  spectral  range  (300  to  900  nm).  This  instrument 
will  be  space  qualified  for  Shuttle  flight.  We  will  seek  a  Shuttle  manifest  slot  in 
collaboration  with  our  subcontractors  at  University  of  Arizona. 

•  Development  of  a  GLO  and  FTHSI-based  instrument  which  will  be  compatible  with  the 
MightySat  II  spacecraft  platform.  We  will  seek  a  MightySat  II  satellite  development  and 
launch  slot  in  collaboration  with  the  USAF  Research  Laboratory  Spacecraft  Division. 

•  Exercise  of  the  developmental  prototype  instruments  on  an  airborne  platform  to 
investigate  their  capabilities  for  detection  of  CAT. 
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6. 


CONCLUSIONS 


We  summarize  the  work  performed  and  the  accomplishments  on  this  SBIR  contract  as 

follows: 

•  Partial  simulations  of  the  background  Gaussian  correlated  random  noise  structure  allowed 
measurement  of  the  PSD  and  Phase  Spectra  and  velocity  for  the  virtual  triangulation 
geometry  concept.  High  signal  to  noise  ratio  in  these  simulations  and  subsequent 
measurements  is  required  to  accomplish  very  high  spatial  resolution  atmospheric  structure 
measurements  from  satellite  altitudes. 

•  More  detailed  simulations,  involving  a  true  3-D  correlated  random  background  drifting 
through  the  triangulation  sensors’  fields  of  view  are  required  and  could  not  be  performed 
in  this  phase  because  of  computer  memory  limitations. 

•  Ground-based  CCD  camera  measurements  made  to  simulate  tropospheric-stratospheric 
structure  measurements  show  the  importance  of  careful  control  of  focal  plane  noise  and 
pattern  noise.  Focal  plane  cooling  will  certainly  be  required  for  visible,  UV,  and  MWIR 
atmospheric  structure  sounding  using  the  techniques  investigated  in  this  project. 

•  The  feasibility  of  using  Fourier  Transform  spectrometers  to  detect  airglow  and  (by 
inference)  atmospheric  backscatter  structure  was  demonstrated  by  the  Kestrel  FTHSI 
instrument  operated  against  an  airglow  background.  Again,  the  importance  of  focal  plane 
and  instrument  cooling  to  minimize  background  noise  was  demonstrated. 

•  Measurements  of  airglow  structure  from  the  Shuttle  using  the  GLO  instruments  shows 
that  the  virtual  triangulation  technique  is  feasible.  Instrument  line-of-sight  ambiguities  and 
spacecraft  attitude  changes  limited  our  ability  to  obtain  detailed  pixel  by  pixel  correlation, 
but  the  overall  data  (integrated  over  the  focal  plane  width)  verifies  that  the  virtual 
triangulation  technique  works. 

•  The  results  of  this  limited  effort  are  promising  but  not  definitive.  Therefore,  more  detailed 
simulations  followed  by  ground-based  experimental  measurements  using  the  airglow 
layers  as  the  signal  source  are  recommended  before  proceeding  to  a  formal  instrument 
design  phase. 
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8. 


APPENDIX  1 


Kestrel  Corporation  Experiment  Report 


This  report  summarizes  the  ground-based  AIRCAM  and  FTHSI  experiment  operating 
parameters  as  conducted  by  the  Kestrel  Corporation. 
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24010078  <  Blue  Camera  Lens  Serial  Number  > 

<  Blue  Channel  V2  > 

412-525,455  <  Yellow  Camera  Filter  412-525, hp-455  > 

002817  <  Yellow  Camera  Serial  Number  > 
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<  Green  Channel  V2  > 
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0010  <  IR  Camera  Serial  Number  > 

22240-084  <  IR  Camera  Lens  Serial  Number  > 
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9. 


APPENDIX  2 


University  of  Arizona  Data 


This  is  the  final  report  from  Dr.  Lyle  Broadfoot,  University  of  Arizona  on  GLO 
experiment  data  analysis. 
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GLO  5  and  GLO  6  Data 
Notes  on  Data  in  following  pages 


Page  1  -  GLO  VIEW  -  see  description  in  “Home  Page” 
‘Vww.glo.ipl.ar  izona.edu/glo” 

Page  2  -  see  ^notes  on  page 

Files  are  named  by  the  MET  time  of  their  recording 

3053530A  is  3/05:35:30  MET 

The  alpha  character  represents  the  CCD  source 


A,  B,  C,  D,  E  are  spectrographs 
F,  G,  H  are  imagers 
I  is  a  tracking  image 

Other  files  are  defined  on  the  page 

Page  3  -  @  b.  1st  -  a  list  of  all  records  from  the  b-CCD,  ie:  b  spectrograph 

Page  5  -  The  Geometric  Parameters  in  the  header  of  each”  *.fit”  files 

Page  7  -  A  plot  of  the  elevation  history  in  data  set  3d0535  and  3e0535 

Data  sets  are  named  by  the  name  of  the  first  file  in  the  data  set 

3d0535  means: 

3  rd  day 
d  flight 
05  hrs 
35  min 


GLO  DATA  file  system: 

Flight  GLO #  DAY  Flight 

STS-53  1  _  _ 

STS-63  2  a 

STS-69  3  b 

STS-74  4  c 

STS-85  5  d 

STS-85  5  e 


Page  8  -  TH  (tangent  height)  of  GLO-5,  GLO-6  in  file  3d_e0535.XLS. 
Page  9  -  Hyperspectral  images  of  02(0,0) 


Page  1 0  -  Spectra  and  Hyperspectral  images 

Page  13  -  Plot  of  geometric  data  showing  top  and  bottom  of  the  slit  (corrected 
Altitude) 

Page  14  -  Hyperspectral  images  of  5  emissions 


File  name 

GLO  6 

60H5  1 

6 

OH(5,l) 

6016300 

■6 

01(6300) 

602  0  0 

6 

02(0,0) 

6015577 

6 

01(5577) 

6Na5892 

6 

Na(5892) 

Page  1 5  -  as  above  for  GLO  5 

50H5_1  5  OH(5,l) 

Page  16  -  Constant  altitude  plots.  The  TH  of  59  km  was  selected  since  that 
elevation  was  common  to  both  instruments. 

Note  the  shift  in  maxima 

Page  17-21 

These  plots  show  the  relative  position  of  the  two  traces  @  59kmTH  for 
both  instruments.  The  top  plot  shows  the  original  position  of  the  data  in 
time.  The  bottom  plot  shows  a  shifted  GLO-6  curve  to  match  the  GLO  5 
curve.  The  time  shift  is  noted. 

They  are  in  order  of  time  delay  between  when  the  signature  showed  in  the 
GLO-6  FOV  pointing  ahead  and  GLO-5  pointing  behind  shuttle  track. 

The  difference  in  the  latitude,  longitude  file  was  6.7  min. 
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OT  66- ST 


O  • 


Stars  =  <  1000.0 >  ftLL3 . STR 


This  is  a  reference  list  of  the  files  in  this  Directory  with  some  comments  that  apply  to  all  data  sets 
The  File  is  dir.xls 
Directory  of  E:'.3d0535 


<DIR>  01-10-99  7:52a. 

<DIR>  01-10-99  7:52a .. 


3053530F  FIT 

17.280 

1/10/99  7:45a 
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63.360 

1/10/99  7:45a 

3053530A.FIT 

3053530C  FIT 

63.360 

1/10/99  7:45a 

3053530C.FIT 

3053530B  FIT 

63.360 

1/10/99  7:45a 

3053530B.FIT 

3053530D  FIT 

63.360 

1/10/99  7:45a 

3053530D.FIT 

3053530E  FIT 

63.360 

1/10/99  7:45a 

3053530E.FfT 

3053556F  FIT 

17.280 

1/10/99  7:45a 

3053556F.EIT 

3053556G  FIT 

17.280 

1/10/99  7:45a 

3053556G.FIT 

3053556A  FIT 

63.360 

1/10/99  7:45a 

3053556A.FIT 

3053556C  FIT 

63.360 

1/10/99  7:45a 

3053556C.FIT 

3060603A  FIT 

63.360 

1/10/99  7:47a 

3060603A.FIT 

3060603C  FIT 

63.360 

1/10/99  7:47a 

3060603C.FIT 

3060603D  FIT 

63.360 

1/10/99  7:47a 

3060603D.FfT 

3060603E  FIT 

63.360 

1/10/99  7:47a 

3C60603E.FIT 

30606 18G  FIT 

17.280 

1/10/99  7:47a 

3060618G.FIT 

306061 8F  FIT 

17,280 

1/10/99  7:47a 

306061 8F.FIT 

3060618A  FIT 

63.360 

1/10/99  7:47a 

306061 8A.FiT 

3060618B  FIT 

63.360 

1/10/99  7:47a 

3060618B.FIT 

306061 8D  FIT 

63.360 

1/10/99  7:47a 

3060618D.FIT 

306061 8E  FIT 

63.360 

1/10/99  7:47a 

3060618E.FIT 

3054306B 

63.360 

1/10/99  7:46a 

3054306B 

File  Directory,  All  Files 

864  file(s) 

43,634,534  bytes 

2  dir(s)  509,362,176  bytes  free 

3O_E0535  XIS 

279,552  01-14-99  7:36a  3d_e0535.xls 

XL  worksheets  for  geometric  data 

A  LST 

7.145 

1/12/99  10:19p 

a.lst 

B  LST 

7,435 

1/12/99  10:19p 

b.lst 

C  LST 

7.087 

1/12/99  1 0:1 9p 

c.lst 

D  LST 

7.377 

1/12/99  10:19p 

d.lst 

E  LST 

7.493 

1/12/99  10:1 9p 

e.lst 

F  LST 

7.145 

1/12/99  1 0:1 9p 

f.lst 

G  LST 

6,333 

1/12/99  1 0: 1 9p 

gist 

H  LST 

144 

1/12/99  1 0:1 9p 

h.lst 

1  LST 

144 

1/12'99  1 0: 1 9p 

List 

Lists  of  data  images  from  each  detector  in  GLO  instrument 

UCAT  MAP 

1,396 

1/8/99  3:02p 

ucaLmap 

Ancillary  data  from  JSC  giving  flight  geometric  information 

GETMET  OUT 

27,464 

1/8/99  1:35p 

GETMET. OUT 

A  file  of  MET  times  and  exposure  times  prepared  to  request  UCAT  data 

MAKEHDF  LST 

23.083 

1/11/99  3:55a 

MAKEHDR.LST 

Header  or  each  spectral  image  in  data  set  usualy  merged  with  UCAT  data 

DIR  TXT 

0 

1/14/99  7:54a 

dir.txt 

File  DIR. 

B  XLS 

29,696 

1/14/99  7:51a 

b.xls 

b.lst  in  XL  format 

O2_0  1  HI  FIT 

46.080 

1/12/99  10:38a 

O2_0  1HY.FIT 

BG0  1HY  FIT 

46.080 

1/12/99  10:38a 

GBO  1.FIT 

BG_REST  FIT 

46.080 

1/12/99  10:38a 

BG_REST.FIT 

O2_0_0H'i  FIT 

46.080 

1/12/99  10:38a 

O2_0_0HY.F!T 

OH6_2HY  FIT 

46.080 

1/12/99  10:38a 

OH6  2HY.FIT 

OH5  1HY  FIT 

46.080 

1/12'99  10:38a 

OH5JHY.FIT 

3G8_3HY  FIT 

46.080 

1/12/99  10:38a 

BG8~3HY.F(T 

OH8__3HY  FIT 

46.080 

1/12'99  10:3ea 

OH8_3HY.FIT 

OI6300HY  FIT 

46.080 

1/12'99  10:38a 

OI6300HY.FIT 

Hyperspectral  Images 
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Volume  in  drive  c.  has  no  label 
Volume  Serial  Number  is  4254-0301 


CSfectocy-  of  E:\2GdU535 


3053530B  FIT 

63,360 

3053556B  FIT 

63.360 

3053567B  FIT 

63.360 

nACagna  nr 
JvoOO  t  £0  f  t  i 

ca.  ~ipri 

3053626B  FIT 

63,360 

3053641 B  FIT 

63.360 

3G53703B-  FIT 

63,360 

3653?11S-rff 

irfi 

OCT,  jtxr 

3053726B  FIT 

63,360 

3053741 B  FIT 

63.360 

3053802B-FIT 

63,360 

■mrini  r\n  f-i-r 

1  11 

r — 1  irn 

3053840B  FIT 

63,360 

3053902B  FIT 

63,360 

3053S09&  FIT 

63,360 

30539240  Fff 

CA-  irn 
Cw.JvU 

3053939B  FIT 

63,360 

3054001 B  FIT 

63.360 

3054011  &  FIT 

63,360 

Ofyr-  rif 

OCT7.X*-/ 

3054028B  FIT 

63,360 

30541 00B  FIT 

63.360 

30544088-  FIT 

63,360 

tatj i  *i r~IT 

jvw  r  £ycr  rrr 

M  TP/L 

30541 37B  FIT 

63,360 

30541 59B  FIT 

63.360 

3054207&  FIT 

63,360 

•JKjZt*CE£ III 

CA-'X Cfl- 
5X77  Jw 

3054237B  FIT 

63,360 

3054258B  FIT 

63.360 

3054306&  FIT 

63,360 

30543210  n T 

3054336B  FIT 

63,360 

3054358B  FIT 

63.360 

3O54405&  FIT 

63,360 

305442OD  FIT 

M-'lgn 

XZZj.J'ZXJ 

3054435B  FIT 

63,360 

3054457B  FIT 

63.360 

3054507B  FIT 

63,360 

3054534B  FIT 

63,360 

3054556B  FIT 

63.360 

3054604&RT 

63,360 

3054C1 0D  Pip 

63.3G0 

3054633B  FIT 

63,360 

3054656B  FIT 

63.360 

3054703B  FIT 

63,360 

•WT  474  nn  r~i~r 

Xvvr  7  ou  rtf 

OJ,vw 

3054733B  FIT 

63.360 

3054756B  FIT 

63J360 

3054862B  FIT 

63,360 

30546T7B  Fff 

63r366- 

3054832B  FIT 

63,360 

3054856B  FIT 

63.360 

3054901 B  FIT 

63,360 

36549160  FiT 

C^i,  ‘VTfl 

3054931B  FIT 

63,360 

3054956B  FIT 

63 .360 

3055003B  FIT 

63,360 

305561 8B  Fff- 

ca.'icti- 

3055030B  FIT 

63,360 

3055056B  FIT 

63.360 

3C5510QB.FIT 

63,360 

1/10/93  7:45a 

30535308. 

1/10/99  7:45a 

3053556B. 

1/10199-  7:45a 

3053557B. 

t/70/99- 7:45a- 

365361 2B. 

1/10/99  7:45a 

3053626B. 

1/10/99  7:45a 

3053641 B. 

1/10/99  7:45a 

3Q53703B. 

1/10/99- 7: 45a- 

365371  IB. 

1/10/99  7:45a 

3053726B. 

1/10/99  7:45a 

3053741 B. 

1/10/99  7:45a 

3053802B. 

rrrwtftr  r:4oa 

36598.1  OB. 

1/10/99  7:45a 

3053840B. 

1/10/99  7:46a 

3053902B. 

1/10/99  7:46a 

3053909B. 

1/10/99- 7:46a- 

3653924B. 

1/10/99  7:46a 

3053939B. 

1/10/99  7:46a 

3054001 B. 

1/10/99  7:46a 

305401  IB. 

1/10/99-7: 46a- 

3654Q26B. 

1/10/99  7:46a 

3054028B. 

1/10/99  7:46a 

30541 00B. 

1/10/99  7:46a 

30541 08B. 

1/10/99- 7:46a- 

3654423B. 

1/10/99  7:46a 

30541 37B. 

1/10/99  7:46a 

30541 59B. 

1/10/99  7:46a 

3Q54207B. 

1/16/99-7: 46a- 

3654222B. 

1/10/99  7:46a 

3054237B. 

1/10/99  7:46a 

3054258B. 

1/10/99  7:46a 

3054306B. 

1/16/99  7:46a- 

3654321 B. 

1/10/99  7:46a 

3054336B. 

1/10/99  7:46a 

3054358B. 

1/10/99  7:46a 

3054405B. 

1/10/00  7:40a 

3654^20B. 

1/10/99  7:46a 

3054435B. 

1/10/99  7:46a 

3054457B. 

1/10/99  7:46a 

3054507B. 

1/16/06  7:46e 

3654522B. 

1/10/99  7:46a 

3054534B. 

1/10/99  7:46a 

3054556B. 

1/10/99  7:46a 

3Q54604B. 

1/16/00 -7:40a 

3054§18B. 

1/10/99  7:46a 

3054633B. 

1/10/99  7:46a 

3054656B. 

1/10/99  7:46a 

3054703B. 

1/16/09  7:46a 

365471 8B. 

1/10/99  7:46a 

3054733B. 

1/10/99  7:46a 

3054756B. 

1/10/99  7:46a 

3054802B. 

t/16/9&-?:46a— 

3654817B. 

1/10/99  7:46a 

3054832B. 

1/10/99  7:46a 

3054856B. 

1/10/99  7:46a 

3054901B. 

1/16/00  7:46a 

365491 6B. 

1/10/99  7:46a 

3054931 B. 

1/10/99  7:46a 

3054956B. 

1/10/99  7:46a 

3055003B. 

1/10/99-7:46a- 

365591 8B. 

1/10/99  7:46a 

3055030B. 

1/10/99  7:46a 

3055056B. 

1/10/99  7:46a. 

30551  fY)B. 

-\r\rr-  *  -«  trr-l  r:T 

30551293  FIT 
30551 56B  FIT 
3055 159B  FIT 


3055229B  FIT 
3055256B  FIT 
3055258B  HT 


3055328B  FIT 
3055356B  FIT 
30553588  FIT 

3055428B  FIT 
3055456B  FIT 
3055500B  FIT 

3055527B  FIT 
3055541B  FIT 
3055603B  FIT 


3055625B  FIT 
3055640B  FIT 
3055702B  FIT 

3055725B  FIT 
305574QB  FIT 
30558018  FIT 

3055824B  FIT 
3055839B  FIT 
3055901 B  FIT 
3655908B  FIT- 
30559238  FIT 
3055938B  FIT 
3060002 B  FIT 
30600100  FIT- 
3060022B  FIT 
3060037B  FIT 
30600598  RT 
Jvvw  i  r  u  r  1 1 
30601 22B  FIT 
30601 36B  FIT 
30601588  FIT 

3060221 B  FIT 
306Q236B  FIT 
3060258B  RT 
3060305B  FIT- 
3060320B  FIT 
3060335B  FIT 
3060357B  FIT 
3060405B-RT 
306041 9B  FIT 
3060434B  FIT 
30604588  FIT 

^ furry* nr n  r-t~r 

3060518B  FIT 
306Q533B.FLT. 
30605568  FIT 

JuWAA/JLr  I  I  i 

306061 8B  FIT 
125  ffle(s)_ 
0  dir(s) 


63r366-  HU 
63.360  1/1 

63.360  1/1 

63.360  1/1 

63:366-  ++ 

63.360  1/1 

63.360  1/1 

63.360  1/1 

63-366-  13+ 

63.360  1/1 

63.360  1/1 

63.360  1/1 

63,360  13+ 

63.360  1/1 

63.360  1/1 

63.360  1/1 
63366-  13+ 

63.360  1/1 

63.360  1/1 

63.360  1/1 

63366-  13+ 

63.360  1/1 

63.360  1/1 

63.360  1/1 

CAJ.JOU  ITT 

63.360  1/1 

63J260  1/1 

63.360  1/1 

m  nrn  4 / « , 
wJ.JUJ  ITT 

63.360  1/1 

63.360  1/1 

63.360  1/1 
63.300  -  13+ 

63.360  1/1 

63.360  1/1 

63.360  1/1 

66366-  13+ 

63,360  1/1 

63J360  1/1 

63,360  1/1 

4  I A 

waXJ  ITT 

63.360  1/1 

63.360  1/1 

63.360  1/1 

m  yn  4/4 
ITT 

63.360  1/1 

63J36Q  1/1 

63.360  1/1 
63,366-  -  13+ 

63.360  1/1 

63.360  1/1 

63.360  1/1 

£3,366-  131 

63.360  1/1 

63.360  VI 

63.360  1/1 

4  1 4 

CC77JWT  ITT 

63,360  1/1 

63J3SO.  1/V 

63,360  1/1 

63:366-  13+ 

63,360  1/1 

7 .920.000  bytes 
508  788  786  hytps  hex* 


1/1Q-99-7046 a- 
1/1099  7:46a 
1/1099  7:46a 
1/1099  7:46a 
++er99- 7:46a- 
1/1099  7:46a 
1/1099  7:46a 
1/1099  7:46a 

4_i.«/1_f~M-'  T-  4f>  — 

IT  7  u  vj  rr^Da 

1/1099  7:46a 
1/1099  7:46a 
1/1039  7:46a 

4  :  i  p,  .ry~,  T.  ar*  - 

I'T  *-  ^Oa 

1/1099  7:46a 
1/1099  7:46a 
1/1099  7:46a 

41/8  ./~V\  T  .  a  m.  _ 

rnursn?  rTOfl 

1/1099  7:46a 
1/1099  7:46a 
1/1099  7:46a 

4xn./vt  ~r .  a  _ 

it i u r jj  r . 
1/1099  7:46a 
1/1099  7:46a 
1/10/99  7:46a 

1/7  i  . 

1/10/99  7:46a 
1/10/99  7:46a 
1/1099  7:46a 

4  /  4  /~l  ./-V'  T.  if 

I  /  !  W  -J  -J  I  . 

1/10/99  7:46a 
1/10/99  7:46a 
1/10/99  7:46a 
1/10399  7:46a- 
1/10/99  7:46a 
1/10/99  7:46a 
1/1099  7:46a 
1/10399  7r4ra— 
1/10/99  7:47a 
1/10/99  7:47a 
1/1099  7:47a 
13+0399  7:47a  - 
1/10/99  7:47a 
1/10/99  7:47a 
1/1099  7:47a 
1/10/9S  7:47a 
1/1099  7:47a 
1/1099  7:47a 
1/1099  7:47a 
+3+9399-7t47a 
1/1099  7:47a 
1/1099  7:47a 
1/1099  7:47a 
V1Q/99-7:4?a- 
1/1099  7:47a 
VI 0/99  7/47a_ 
1/1099  7:47a 
W+999-?r4?a- 
1/1099  7:47a 
1/1899  7:47a 
1/1099  7:47a 
13+9399-7:4rs- 
1/1099  7:47a 


36591 15B. 
3055 129B. 
30551 56B. 
3055159B. 
3655?14B. 
3055229B. 
3055256B. 
3055258B. 
365531 3B. 
3055328B. 
3055356B. 
3055358B. 
365541 3B. 
3055428B. 
3055456B. 
3055500B. 
365551 4B. 
3055527B. 
3055541 B. 
3055603B. 
3655§11B. 
3055625B. 
3055640B. 
3O55702B. 
36S5710B. 
3055725B. 
3Q55740B. 
3055801 B. 
3655809B. 
3055824B. 
3055839B. 
3055901 B. 
3655§08B. 
3055923B. 
3055938B. 
3Q60002B. 
366001 OB. 
3060022B. 
3060037B. 
3060059B. 
366Q107B. 
30601 22B. 
30601 36B. 
3060158B, 
366Q206B. 
3060221B. 
3060236B. 
3060258B. 
366Q305B. 
3060320B. 
3D60335B. 
3Q60357B. 
306Q405B. 
3060419B. 
3060434B. 
3060458B. 

n/wvTnrn 

v/WA/wfUQu. 

306051 8B. 
3DS0533B. 
3060556B. 
366Q603B. 
306061 8B. 
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3e0535  01(6300) 


■GLO  S  53 fcm 


’GLO'&4-2.0fnffl 


Minutes  from  3/05:30 


3e0535 


Geometric  Data;  Set  3d0535 


Geometric  Data;  Set  3e0535 
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c imos/WM 


Vanguard  Research,  Inc. 


www.vriffx.com 

10400  Eaton  Place,  Suite  450 
Fairfax,  VA  22030-2201 
Phone  (703)  934-6300 
Fax  (703)  273-9398 


May  4,  2000 
00-0292 


Defense  Technical  Information  Center 
8725  John  J.  Kingman  Road 
Suite  0944 

Ft.  Belvoir,  VA  22060-6218 

Subject:  Final  Report  -  Replacement  Page,  SF298 

Reference:  Contract  No.  DASG60-98-M-0073 


Dear  Sir  or  Madam: 

Please  find  enclosed  two  copies  of  a  replacement  page  for  SF298  for  the  final  report.  Please  replace 
this  page  in  your  copies  of  the  final  report. 


Should  you  have  any  questions,  do  not  hesitate  to  contact  me  at  (703)  934-6300. 


Sincerely, 


Debra  A.  Spear 
Contracts  Administrator 


Enclosure  as  stated 


